Electron Paramagnetic Resonance of Boron Acceptors in Isotopically Purified Silicon 
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The electron paramagnetic resonance (EPR) linewidths of B acceptors in Si are found to reduce 
dramatically in isotopically purified ^*Si single crystals. Moreover, extremely narrow substructures 
in the EPR spectra are visible corresponding to either an enhancement or a reduction of the ab- 
sorbed microwave on resonance. The origin of the substructures is attributed to a combination of 
simultaneous double excitation and spin relaxation in the four level spin system of the acceptors. A 
spin population model is developed which quahtatively describes the experimental results. 



Boron is the most widely employed shallow acceptor in 
silicon. However, due to the existence of four degenerate 
energy levels originating from the light and heavy hole 
bands (shown in Fig. [l](a)) and their strong sensitivity 
to perturbations like internal strain [IJ, B exhibits rich 
optical and electron paramagnetic resonance (EPR) spec- 
tra, which are still not fully understood. Recently, some 
open questions associated with the optical properties of 
this acceptor have been solved due to the availability of 
isotopically enriched Si 0,11]. Using such samples, it has 
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FIG. 1: (a) The level schemes of B acceptors. An externally 
applied magnetic field splits the heavy (light) hole subband 
into two energy states |±3/2) (|±l/2)). We refer to each 
level as |1)-|4) in the order shown. On top of the EPR-allowed 
transitions defined by Am=l, the experimentally observed 
nominally EPR-forbidden transitions Am=2 and Am=3 are 
also shown, (b) EPR spectra recorded at r=3.0 K on ^*Si:B 
and "^'Si:B samples showing the Am=l and 2 transitions. 



been shown that the random distribution of the different 
host Si isotopes causes a residual acceptor ground state 
splitting in natural Si which is experimentally observed 
in the photoluminescence (PL) spectra of acceptor bound 
excitons 4] and also gives rise to inhomogeneous broad- 
ening of the infrared absorption spectra 0, 0| • 

Similarly, host isotope effects might also account for 
various properties of the EPR spectra of B in Si that 
have not been understood so far. The first successful 
boron EPR required the application of large homoge- 
neous external stress to diminish the influence of inho- 
mogeneous internal strain flelds that are mainly induced 
by randomly distributed point defects and dislocations. 
These defects lead to a statistical distribution of ground 
state splittings, which, depending on the defect concen- 
tration, can result in an extremely large inhomogeneous 
broadening of the EPR resonance lines in externally un- 
strained Si [7| . The first observation of boron EPR with- 
out external stress was finally realized in 1978 , when 
Si crystal growth techniques had developed enough to 
obtain samples of sufficient crystalline quality. However, 
the experimentally observed EPR linewidths have always 
been much broader than could be explained solely by 
the residual defect concentrations in high quality Si. In 
addition, Neubrand observed sharp substructures in the 
main resonance peaks, with and without externally ap- 
plied stress d, . More recent studies on acceptor EPR 
probing dependencies on temperature, external stress, 
and frequency [3, 11 1 found the same sharp features. 



but the origin of these substructures remained unclear. 

Here, we present experimental evidence that indeed 
the random distribution of stable Si isotopes is the ma- 
jor cause of inhomogeneous broadening of boron EPR 
spectra in natural Si. Moreover, thanks to the sharpen- 
ing of the main EPR lines, the improved visibility of the 
sharp substructures allows for detailed studies of their 
origin. The substructures are shown to arise from an 
interplay between simultaneous excitation of two EPR 
transitions and spin relaxation in the four level spin sys- 
tem of B acceptors in Si. A theoretical model based on 
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this assumption is developed to describe the temperature 
dependence of the substructures. 

The experiments were performed with two p-type Si 
single crystals. The first sample, referred to as "'^^SiiB, 
was float-zone grown, has a natural isotopic composi- 
tion and a B concentration of 1.4x10^'' cm~'^. The iso- 
topically purified sample (~99.98 % ^*Si), referred to as 
^*Si:B, comes from the neck of a float-zone crystal and is 
doped with B to a concentration of 3(±l)xl0^'* cm~'^. 

EPR measurements at a temperature of T=1.8 K 
were performed using a Bruker E580 EPR spectrome- 
ter equipped with a flex-line resonator (ED-4118MD5) 
operating at a microwave frequency i>=9.6 GHz. An Ox- 
ford CF935 helium-flow cryostat was used to maintain 
the temperature. For experiments at T>2.8 K, a Bruker 
E500 spectrometer in conjunction with a super-high-Q 
resonator (ER-4122SHQE) was used at i/=9.4 GHz. In 
this case, the sample was cooled using an Oxford ESR900 
helium-flow cryostat. Lock-in detection using magnetic 
field modulation leads to the detection of the first deriva- 
tive dx" /dH of the imaginary part of the magnetic sus- 
ceptibility x" with respect to the magnetic field H . The 
absorption curves, which are proportional to x"j are ob- 
tained via numerical integration of the measured spectra. 

Figure [1] (b) shows two overview EPR spectra that 
were measured on the ^^Si:B and the "'^^SiiB samples. 
In both spectra, we observe three distinct B-related res- 
onances. From the angular dependence of the effective 
g- values [1,01, we assign the resonance at iJ=2980 G to 
a superposition of the two inter-subband Ato=2 tran- 
sitions, where m is the magnetic quantum number of 
the total angular momentum J=3/2 of the B acceptor. 
The resonances at iJ=5730 G and i?=6220 G originate 
from a superposition of the two inter-subband lS.m—1 
transitions |1)-(-t>|2) and |3)-f->-|4) and the intra-subband 
Am=l transition |2)-(^|3), respectively [1, i, The 
most prominent difference between the two spectra is 
the decreased linewidth of all B-related resonances in the 
isotopically purified sample. This effect is most strongly 
pronounced for inter-subband transitions. As found by 
PL in Ref. 0, @] , local fluctuations of the valence band 
edge due to the random distribution of the different Si 
isotopes are responsible for the residual splitting of the 
B acceptor ground state in "'^'Si in absence of an ex- 
ternal magnetic field. A detailed modelling of the in- 
fluence of the isotope-induced fluctuations on the EPR 
resonances shows that this effect also accounts quanti- 
tatively for the inhomogeneous broadening of B-related 
EPR signals that could not be explained by point defect- 
induced random strain fields. This is in marked differ- 
ence to the broadening of donor EPR lines in Si, which 
is caused by superhyperfine interaction with ^^Si nuclei 
[l3| . The quantitative investigation of the broadening 
mechanism in Si:B is beyond the scope of this paper and 
will be given elsewhere. Here, we focus on the origin of 
the sharp substructures that appear in the middle of the 



inter-subband resonances. The observation of these sub- 
structures in ^^Si:B establishes that they do not result 
from an isotopic effect. 

Figure [2] shows the inter-subband resonances on an ex- 
panded magnetic field scale as measured for ^^Si:B at 
two different temperatures. The narrow substructure 
is observed in all four spectra. By cooling from 3.2 to 
1.8 K, the linewidth of the substructures decreases sub- 
stantially, from 15 to 3 G for Am=l and from 5 to 1 G 
for Am=2. In contrast, the broad underlying signals 
show much weaker changes in linewidth with tempera- 
ture (from 78 to 70 G for Am=l, 36 to 35 G for Ato=2). 
Thus, even though the line shape of the broad resonances 
is Lorentzian, the lack of any significant temperature de- 
pendence suggests that the lines are inhomogeneously 
broadened. This broadening is attributed to the dis- 
tribution of transition energies between different energy 
levels induced by random local strains, which are well ac- 
counted for by the ~1 x 10^^ cm~'^ concentration of C im- 
purities as measured by infrared absorption spectroscopy 
on the ~607 cm~^ C local vibrational mode [jl. [l^. 

The insets in Fig. [2] show the integrated absorption 
curves in the magnetic field region of the substructure. 
We observe the following experimental facts: (1) For the 
Am=l spectra (Fig. [2] (a)), the substructure has a neg- 
ative sign corresponding to a decrease of the microwave 
absorption at the center of the broad resonance. The in- 
tensity of this dip increases with decreasing temperature. 
We have exclusively observed a negative sign of the sub- 
structure in the Am=l resonance for all temperatures 
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FIG. 2; EPR spectra of the inter-subband Am=l resonance 
(a) and the Am=2 resonance (b) in ^*Si:B for -WHflOO]. The 
upper spectrum in each figure was measured at T=3.2 K, 
v—Q.SS GHz and a microwave power of P=0.2 mW. The 
lower spectra were measured at r=1.8K, u—Q.SG GHz and 
P=0.002 mW. The values of P were chosen to achieve the 
same microwave excitation rate w in the two resonators used. 
Solid lines show the results of numerical fittings composed of 
two Lorentzian lines, which are separately shown as dashed 
and dotted lines. Ho is the resonance field of each signal. The 
amplitudes of the spectra were normalized. The insets show 
the integrated absorption spectra in the magnetic field region 
of the substructures in arbitrary units. 
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investigated. (2) In contrast, for the Am=2 resonance, 
the microwave absorption in the substructure is enhanced 
at high temperatures, giving rise to a sharpened peak of 
the integrated spectrum measured at T=3.2 K. However, 
when the temperature is decreased to 1.8 K, the sign of 
the substructure hue reverses, and a dip in the absorp- 
tion spectrum is observed. (3) Microwave power satura- 
tion experiments (not shown) reveal that the saturation 
behavior of the substructure lines strongly deviates from 
that of the underlying broad signals; the relative intensity 
of the substructure decreases with increasing microwave 
power. No power broadening is observed for the substruc- 
ture lines. (4) As has been pointed out in Ref. 0, the res- 
onance fields of the substructure lines do not change when 
a small external stress is applied. (5) We find that the 
linewidths of the intra-subband Ar7i=l resonance and the 
substructure of the inter-subband Am—1 signal are very 
similar, in particular for temperatures below 3 K. More- 
over, in the "^*Si sample, for each orientation of H (not 
shown) these two resonances have the same asymmetry 
of the line shape, defined as the a/b ratio in Fig. 1(b), in- 
dicating that transitions between |2) and |3) are involved 
at least in the generation of the Ato=1 substructure. 

The experimental observations summarized above 
strongly argue against a model assuming that the nar- 
row substructure lines arise from a partial overlap of two 
asymmetric lineshapes, e.g. the |l)o|2) and |3)^|4) 
transitions in case of the Ato=1 signal, or the |l)o|3) 
and |2)o|4) transitions in case of the Am=2 signal, as 
has been already speculated in Ref. [1]. On the other 
hand, a similar "substructure" line in the middle of a 
broad resonance has been reported for Ni^+ centers in 
MgO [15]. However, the cross-relaxation model invoked 
there can not account for the positive sign of the Am=2 
substructure observed here for T>3.2 K. In what follows, 
we demonstrate which interplay of spin excitation and re- 
laxation in the four level system of Si:B causes the narrow 
substructure lines observed. 

The centers of the broad resonances, where the sub- 
structure lines develop, correspond to B acceptors that 
experience a strain close to zero. In this case, the en- 
ergy level scheme of the four level system is symmetric, 
as shown in the insets of Fig. |3l and the energies of the 
transitions |l)o|2) and |3)o|4) [|l)o|3) and |2)o|4)] 
are equal, so that they are resonantly excited at the same 
time. This double excitation can lead to a change in the 
microwave absorption with respect to the situation when 
only one transition can be excited. The magnitude and 
sign of this change will depend on the ratios between 
relaxation and excitation rates in the level system. 

To explore this scenario, we have performed a first 
calculation based on an incoherent spin population 
model composed of a set of rate equations. The 

time evolution of the spin population n, of each state \i) 
is given by the sum of all possible excitation and relax- 
ation terms. Relaxation and excitation rates from \i) to 
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FIG. 3: Contour plots obtained by the population model de- 
scribed in the text for (a) Am—1 and (b) Am=2. The contour 
lines indicate a relative intensity (in %) of the substructure 
lines as compared to the intensity of the broad underlying 
transition. Negative and positive intensity corresponds to 
a suppression (dip) and an enhancement (peak) of the mi- 
crowave power absorptions, respectively. The negative region 
is shown shaded. The parameters used in the calculation were 
r=3 K, p=9A GHz, and 714/^=712/^=1, with w being the 
resonant microwave excitation rate. The dashed line repre- 
sents the assumed correlated temperature dependence of 713 
and 723 in ^*Si:B as discussed in the text. 



|j) are represented by 7^^ and Wij, respectively. 

We use our model to compare the microwave absorp- 
tion between exactly on and slightly off resonance con- 
ditions. For the on resonance condition where the sub- 
structure emerges, W12 and W34 [wis and W24] are as- 
sumed to have the same non-zero value w for the Am=l 
[Ar7i=2] resonance, while all the other elements in the 
excitation matrix are set to zero. The rate equations are 
solved for the steady state condition (dn/dt^O) using the 
following assumptions: (i) Boltzmann balance of the re- 
laxation rates, i.e. "1ji=lij exp {—AEij /kT), where AEij 
is the energy difference between the states \i) and |j) and 
k is the Boltzmann constant, (ii) 712=734 and 713=724, 
dictated by the symmetry of the Si:B Hamiltonian when 



the local strain is close to zero, (iii) Wji 



required for 



magnetic dipole transitions, and (iv) Ui—l to conserve 
the total number of spins. The sum of the resulting popu- 
lation differences, {ni-n2)+{n3-n4) [{ni-n3)+{n2-n4)] for 
the Am=l [for Ato=2] transition, is used as a measure 
for the microwave absorption. To simulate the situation 
slightly off resonance, where double excitation is not pos- 
sible, only one excitation at a time is allowed, e.g. Wi27^0 
[wi3 7^0] is assumed for the Am=l [Am=2] resonance, 
and the microwave absorption represented by ni-n2 [ni- 
n^l is added to the microwave absorption n^-n^ [n2-n4] 
that is obtained by separately considering only the sec- 
ond transition with the excitation w^i^iO [w247^0]. 

Results of such calculations are shown in contour plots 
in Fig. [3) The relative change of the microwave absorp- 
tion due to double excitation, calculated as the popula- 
tion differences between exactly on and slightly off reso- 
nance conditions, is plotted as a function of the normal- 
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ized relaxation rates 713/u; and 723/^1 (all other relax- 
ation rates were kept constant, e.g. ji2/w=ju/w=l). 
This calculated change in absorption can be compared 
with the intensity of the substructure in Fig. [5] The open 
and shaded regions in Fig. [3] correspond to a suppression 
(dip substructure) and an enhancement (peak substruc- 
ture) of microwave absorption, respectively. The dashed 
line shown in the plot corresponds to a constant ratio of 
713/723=4/1, motivated by the temperature independent 
ratio of substructure linewidths, Ai7Am=i/Ai?Am=2~4, 
observed experimentally. Assuming all other parameters 
to only depend moderately on temperature, a variation 
of the temperature in our experiments can be viewed in a 
first approximation as a correlated change of 713 and 723 
along the dashed line. This picture qualitatively explains 
the experimentally observed temperature dependence of 
the substructure line intensity in Fig.[2l For example, we 
find a negative sign (dip) for the Am=l substructure and 
a positive sign (peak) for the Am=2 substructure when 
7i3/w>4 and 723/w>l. Experimentally, this situation 
corresponds to temperatures >3.2 K. Upon decreasing 
both 713/w and 723 /w along the dashed line, correspond- 
ing to lowering the temperature from 3.2 K to 1.8 K in 
the experiment, we observe a change of the sign for the 
Am=2 substructure from positive to negative. For the 
Am=l substructure, the sign always stays negative for 
all 7i3 and 723 along the dashed curve, agreeing with the 
experimental finding that the sign is always negative at 
all accessible temperatures. 

By simulating the model using a range of jij param- 
eters we have identified several requirements for the rel- 
ative magnitude of the 7.^ in order to reproduce the ex- 
perimental observations. First, to achieve a substantial 
effect (e.g. 10 %) from the double excitation requires that 
all 7ij are of the same order of magnitude and also that 
the applied microwave power is close to saturation, i.e. 
7i2~7i3'^7i4~723'^'w, as used in Fig. [3l Second, to ob- 
serve a temperature dependent change of the sign for the 
Am=2 substructure requires that 7i3>7i2. This order of 
the jij is unusual and contrasts with 7i2>723S>7i3S>7i4 
usually seen in other high-spin systems where zero-field 
splitting (ZFS) dominates spin relaxation However, 
we point out that the ZFS Hamiltonian in Si:B is rep- 
resented by the Sf terms [8j] (in contrast to usual Sf ), 
and therefore the order of relaxation rates can be differ- 
ent. The results presented here should therefore moti- 
vate the development of a detailed theory of relaxation 
effects in Si:B. Independent of this, we have performed 
microwave power saturation measurements (not shown) 
and confirmed that all resonances in Si:B saturate at 
about the same microwave power. The multiple relax- 
ation pathways in the four level system make a direct 
interpretation of these saturation experiments difficult. 
Nevertheless, the fact that microwave power saturation 
occurs at the same power indicates that at least some of 
the jij in Si:B have comparable magnitudes as assumed. 



Our model is still somewhat qualitative and does not 
reproduce all experimental facts. For example, the model 
predicts that the strongest effect from double excita- 
tion should be seen at close to saturation conditions, i.e. 
"/ij/w^l, which is not in line with experiment. Further 
improvements to the model could be based on the density 
matrix formalism and the Redfield relaxation theory to 
construct the full set of Bloch equations for the J=3/2 
spin system, also taking into account the specific spin 
Hamiltonian of Si:B to correctly predict the relaxation 
rates between difi'erent spin levels. Nevertheless, the in- 
structive model constructed here establishes that the sub- 
structures in the EPR lines arise from the double excita- 
tion among the four levels of B acceptors. The improved 
understanding of the resonance properties of acceptors 
and the additional degree of freedom achieved via iso- 
tope engineering could make these valence band-derived 
states versatile competitors to the more thoroughly stud- 
ied donors also in applications such as spin- or orbitron- 
ics. 
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